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We consider the exact reduced dynamics of a two-level system coupled to a bosonic reservoir, 
further obtaining the exact time-convolutionless and Nakajima-Zwanzig non-Markovian equations 
of motion. The considered system includes the damped and undamped Jaynes-Cummings model. 
The result is obtained by exploiting an expression of quantum maps in terms of matrices and a simple 
relation between the time evolution map and time-convolutionless generator as well as Nakajima- 
Zwanzig memory kernel. This non-perturbative treatment shows that each operator contribution 
in Lindblad form appearing in the exact time-convolutionless master equation is multiplied by a 
different time dependent function. Similarly, in the Nakajima-Zwanzig master equation each such 
contribution is convoluted with a different memory kernel. It appears that depending on the state of 
the environment the operator structures of the two set of equations of motion can exhibit important 
differences. 
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I. INTRODUCTION 



The study of open quantum systems is a wide research 
area of interest to various scientific communities, rang- 
ing from physicists to chemists and mathematicians. Its 
basic theoretical framework is well-understood when a 
Markovian description can be applied [l|, [2| , but despite 
important work [2j, |3j a lot remains to be clarified and 
understood when considering non-Markovian dynamics, 
which is often mandatory for a realistic approach, also 
coping with strong coupling between system and envi- 
ronment. Recently a lot of research work has been de- 
voted to the subject |4l-[24j. On the one hand it is not 
fully clear which is the most general operator structure 
of non-Markovian equations of motion, which do provide 
a well-defined time evolution and in particular preserves 
complete positivity. On the other hand one would like 
to link, in a possibly intuitive way, operator structures 
giving a sensible dynamical evolution with microscopic 
information on the physics of the system of interest. For 
the Markovian case both these approaches, phenomeno- 
logical and microscopic, are well understood and success- 
fully used. The result by Gorini, Kossakowski, Sudarshan 
and Lindblad [25l l26j provides a robust framework in or- 
der to envisage a possible ansatz for the evolution equa- 
tions, while the Markov approximation is often connected 
with weak coupling, thus making microscopic approaches 
more manageable. In considering non-Markovian dynam- 
ics, there is a natural tendency to bring over the intu- 
ition gained for the Markovian case. As it appears, how- 
ever, the non-Markovian case is much more subtle and 
involved, so that in such a way one often runs into incon- 
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sistencies and pitfalls. 

In this article we will discuss a realistic physical model, 
simple enough to be exactly treated in detail, but already 
showing typical non-Markovian features. In particular 
we not only derive the exact equations of motion, which 
provide the counterpart to the Lindblad structure of a 
Markovian model, but also obtain a clear-cut connec- 
tion to the microscopic approach via the two standard 
Nakajima-Zwanzig and time-convolutionless techniques. 
This work should help in better understanding typical 
non-Markovian features in simple but realistic models. 
We will consider a two-level system coupled first to a 
single mode of the radiation field, and later to a bath of 
harmonic oscillators, via a Jaynes-Cummings type of in- 
teraction. By exploiting the knowledge of the exact uni- 
tary evolution, and therefore of the reduced dynamics, as 
well as a suitable matrix representation of the dynami- 
cal maps, we can exhibit the exact time-convolutionless 
and Nakajima-Zwanzig equations of motion. It will ap- 
pear that such equations are given by a sum of Lindblad 
terms, each multiplied or convoluted with a different time 
dependent function. Moreover the operator structure of 
these two sets of equations will be shown to be gener- 
ally different, also depending on the environmental state. 
These exact results, with a clear-cut connection to a mi- 
croscopic perturbative derivation, should guide in devel- 
oping a correct intuition for the approximate treatment 
of more involved non-Markovian systems. 

The article is organized as follows. In Sect. |IT]we de- 
rive the map giving the exact reduced dynamics for the 
Jaynes-Cummings model, further introducing a simple 
representation of the map as a matrix with respect to 
a convenient basis of operators. In Sect. IIIII we point 
out the simple relations connecting the time evolution 
map to the time-convolutionless, time-dependent gener- 
ator and the Nakajima-Zwanzig memory kernel, respec- 
tively. These relationships are most easily expressed in 



matrix form. The corresponding matrices are evaluated 
for the reduced dynamics of our two-level system, fur- 
ther coming back to the operator expression of the mas- 
ter equations and giving the explicit result for a bath in 
the vacuum state. In Sect. IIVI we put forward a similar 
analysis for the damped two-level system, making contact 
with previous results and comparing with the undamped 
case. 



II. JAYNES-CUMMINGS MODEL AND EXACT 
REDUCED DYNAMICS 

A. Time evolution mapping 

We consider a two-level system coupled to a single 
mode of the radiation field according to the total Hamil- 
tonian 



where the following operators have been introduced 



H — Hs + He + Hj, 

where the system Hamiltonian is given by 
H s = ui a+cr-, 



(1) 



(2) 



with wo the transition frequency, a+ = |1)(0| and er_ = 
|0}(1| the raising and lowering operators of the two- level 
system. The Hamiltonian for the single mode of the ra- 
diation field is given by 



H h 



U! 



tfb. 



(3) 



where the creation and annihilation operators b* and b 
obey the standard bosonic commutation relation. The 
coupling is in the Jaynes-Cummings form 



Hr 



g(a + <E)b + a-<g) b*) , 



(4) 



so that the considered model can describe e.g. the in- 
teraction between a two-level atom and a mode of the 
radiation field in electric dipole and rotating wave ap- 
proximation. Working in the interaction picture with 
respect to the free Hamiltonian Hs + He, 



with 



Hi(t) = g (cr+ ® be lAt + ct_ ® &V iAt ) , (5) 



A = LUq — U! (6) 



the detuning between the system and the field mode, it 
is possible to obtain the exact dynamics generated by 
the total Hamiltonian (see e.g. Puri |27|), and there- 
fore the reduced dynamics of the two-level system. We 
express the result exhibiting the unitary evolution oper- 
ator which in the basis {|1}, |0}} is given by the following 
matrix, whose entries are operators in the Fock space of 
the radiation field 



U(t) = 



c(h+l,t) d(h+l,t)b 
-6 t d + (fi + l,t) c" (n,t) 



(7) 



c{n,t) 



JAt/2 



cos I \/A 2 + 4g 2 n- ) 

sin [\/A 2 + 4g 2 nl- 

sA V ; — 

V A 2 + Ag 2 h 



(8) 



d (n, t) 



-le 



iAt/2 



2ff- 



sin ( \/A 2 + 4g 2 n| 



a/A 2 + 4g 2 h 



(9) 



with h = b*b the number operator. The unitarity of U (t) 
is granted because of the easily verified relation 



c* (n,t) c(h,t) + nd* (n,t) d(n,i) = 1. 



(10) 



Given the unitary evolution of the whole bipartite system 
one can obtain the reduced dynamics of the two-level 
atom simply by taking the partial trace with respect to 
the environmental degrees of freedom. If the initial state 
of system and environment is factorized, the map giving 
the reduced dynamics is completely positive and takes 
the form 



p(t) = Tr B [U(t)p(Q) ®p E U\t)] 



(ii) 



Taking U(t) as in Eq. (J7J) and considering an envi- 
ronmental state commuting with the number operator, 
[p£;,n] = 0, so that in particular both the vacuum and a 
thermal state can be dealt with, one comes to the follow- 
ing explicit expression for the action of the map $(£): 

Pii(t) = Poo (0)(l-a(t))+pn (0)0(4), 

Pio(t) = Pio(0) 7 (i), 

Poi(t) = Poi(0) 7 *W, 

Poo(i) = poo (0) a (t) + pu(0)(l -/?(£))■ (12) 

The effect of the interaction with the bath is contained 
in the time dependent coefficients a(t), j3(t) and "f(t), 
which are given by the following expectation values over 
the state of the environment pe '■ 

a(t) = (c f (h,t)c(h,t)) E , 

f3{t) = (c t (n + l,t)c(n + l,t)) Bl (13) 

7 (t) = (c(n,t)c(h + l,t)) E - 

B. Matrix representation 

Now that we have obtained the completely positive 
map <&(i) giving the exact reduced time evolution of the 
considered two-level system, we shall express it for later 
convenience in matrix form with respect to a suitable 
basis of operators in the Hilbert space (D 2 , following [28| ■ 
Indeed, having fixed a basis {Xi} l=0 1 2 3 of operators in 



(D 2 , orthonormal according to Tr< 

map A can be expressed in this basis according to 



*l*l 



Ski , a linear 



Ap = J2 L ki^s\xjp 



Xh 



(14) 



1,1 



where the matrix of coefficients Lhi uniquely associated 
to the map is given by 



Lki = Trc 



xlHXi) 



(15) 



Here and in the following we will use Greek or calligraphic 
letters to denote maps, and Roman letters to indicate the 
corresponding matrix. The convenient choice of basis for 
our calculations, in order to later recast the relevant maps 

in operator form, is given by {Xi} l=Qi = |^1, -^ a i\> 
where cr., denote the usual Pauli operators. This choice 
leads to the following expression for the matrix Fki(t) 
associated to the time evolution map $(£): 



F(t) 



m 



1 















\ 







7fl(*) 


li (t) 














-li (t) 


m(t) 







' 


— 


a(t) 








Pit) 


+ a(t) 


-1/ 

(16) 


1C1 


ents 


defined 


in Eq. (TT5} i 


ippear, 


and we 



denote with R and / real and imaginary part of a given 
function: 

1 = 1r + ili- 

This strategy of associating matrices to maps, already 
pursued in .28] , and which transforms composition in ma- 
trix multiplication, will turn out to be very convenient to 
obtain the expressions of exact equations of motion lead- 
ing to the dynamics described by Eq. (fi"2"|) . 



III. EXACT NAKAJIMA-ZWANZIG AND 

TIME-CONVOLUTIONLESS MASTER 

EQUATIONS 

A. General expression in terms of time evolution 
map 

With the aid of the knowledge of the exact time evo- 
lution, and using the representation of maps in terms of 



matrices, we will now explicitly obtain two kinds of exact 
equations of motion for the reduced system's dynamics. 
We first consider a master equation in differential 
form with a generator local in time, i.e. the time- 
convolutionless master equation. Assuming the existence 
of such a generator /Ctcl (t) , it should obey the equation 

P(t) = KtclCOpC*), ( 17 ) 

which due to p(t) = $(£)/? (0) is satisfied upon identifying 

Ht)§- l {t), (18) 



£tcl(*) 

or in terms of matrices 



^tclW = F(t)F- x (t), 



(19) 



holding provided the inverse does exist. 

On a similar footing one can consider a master equa- 
tion in integrodifierential form, given by a suitable mem- 
ory kernel, corresponding to the Nakajima-Zwanzig mas- 
ter equation. In this case the memory kernel /Cnz(^) 
should obey the convolution equation 

p{t) = (£nz *P)(*), (20) 

so that in view of Eq. (jlll) one has the relation 



£nz(u) = ut-Q-^u), 



(21) 



where the hat denotes the Laplace transform, and there- 
fore in matrix representation 



K NZ {u) = ul-F-^u). 



(22) 



It immediately appears that, given the time evolution 
map from the relationships Eq. (fH))) and Eq. (|2"2"j) , one can 
directly obtain the generator of the master equation in 
time-convolutionless form, or the memory kernel for the 
Nakajima-Zwanzig form respectively, without resorting 
to the evaluation of the whole perturbative series, which 
of course leads to the same result [2(| . Obviously, given 
the exact time evolution one does not need the equations 
of motion. Nevertheless an exact comprehensive study as 
feasible in this case allows, as we shall see, to point out 
general features, serving as a guide for phenomenological 
or approximate treatments. 

Starting from Eq. (|16|) and Eq. (fT9|) one obtains for the 
model of interest 



J 



!<■ 



TCL 



(*) 









[l-2/3(t)]q(t)-[l-2a(t)]/3(t) 
P(t)+a{t)-l 



Re 
-Im 





" id 

7(t) 

7(*) 









Im 
Re 



7(*) 
7(t) 
7(0 
7(0 












f3(t)+ a (t)-i 



(23) 



and the expression is well defined provided the determinant 



det F(t) = \j(t)\ 2 [a(t) + /3{t) - 1} 



(24) 



is different from zero. 

On a similar footing one can consider the Laplace transform of Eq. (|16p . given by the matrix F(u) with determinant 



det F(u) 



1R (m)+7/ 2 (m) 



a{u) + j3(u) 



1 



(25) 



and using Eq. 



one further obtains 

/ 
Knz(u) = 



^ u 2 [&(u)-l3(u)} 
\ l-u[a(u)+p(u)\ 



0_ 

_ 7fl(") 

7i(") 
1R 2 (uj+jp (u) 



J) 

7i(«) 

— 7-bC") 

7n 2 («)+77 a (") 







2u-« 2 [q(«)+/3(u)] 
l-u|d(u)+/3(u)| 



(26) 



r 



which upon inverse Laplace transform provides the exact 
Nakajima-Zwanzig integral kernel. As it appears, work- 
ing with the matrix representation has proved very con- 
venient to easily obtain the maps fixing the differential 
and intcgrodifferential equations of motion for the model, 
given by Eq. (fTT|) and Eq. (|2U)) respectively, in terms of 
the completely positive time evolution map Eq. (1121) . 



whose last term can be written in alternative ways ac- 
cording to the identities 



azpa* 



P = 4 

= 4 



a+a-pa+a- - - {a+(T-,p} 
cr-a+pa-a + - - {a^a + ,p} 



(29) 



B. Operator expression 

We now recast the obtained maps, providing time local 
generator and memory kernel, in operator form, to better 
compare with previous work and appreciate the difference 
in the obtained expressions. This is easily done observing 
that a matrix of the form 



Exploiting this result one obtains the exact time- 
convolutionless master equation describing the reduced 
dynamics of a two-level atom coupled according to the 
Jaynes-Cummings model to a single mode of the radia- 
tion field, which is of the form Eq. (fT7|) with /Ctcl (t) 
given by 



.4 





E r Ei 

-Ei E r 

X Y 



(27) 



in the basis < -75 1, "75 °j f corresponds in operator fori 
to the map 



Ap = iEi [a + a-,p] 



- 2 {X + Y) 



(28) 



a +P°- - 2 {v-V+'P} 
< J -P (7 + - 2 {o+0"->P} 



(Y - 2E r ) [a z pa z - p] , 



/Ctcl(*)p = ilm 



7(*) 



Mt) 

[a(t) - 1] $ (t) - /3(t)A(t) 
0(t) + a(t) - 1 

[p(t) - 1] a (t) - a(t)$(t) 
j3(t) + a(t) - 1 

1 \ 0(t) + a{t) 



[a + a^,p] 



°+po- - 2 { a - a +'P} 



4 I pit) + a (t) - 1 



2 Re 



a-pa+ 
7(*) 



-{a+a-,p} 



lit) 



[<*zPO z - P\ 



(30) 



In a similar way one has for the Laplace transform of the 
memory kernel /Cnz(^), appearing in the exact Nakajima- 



Zwanzig master equation Eq. (l20l) . the expression 



the explicit expression 



1Cnz(u)p = i 



II ( u ) r l 

— t. 7. \a + a-,p\ 

■^—2/ \ . ^2/ \ L ~ - rJ 

1R (u) + 7/ ( u ) 



1 



u [ua{u) — 1] 



1-u 



a(u)+ fi{u) 
uj3{u) - 1 



°+P°- - 2 W-V+iP} 



1-u 



a(u)+/3{u) 
&(u)+ $(u) 



a -pa j, 



W+v-,p} 



1r(u) 



a(u) + (3(u) 



1R 2 ( U ) +li 2 {u) 
x WzPcr z - p] . (31) 



Despite being exact these expressions are quite cumber- 
some, since the functions given in Eq. (fT3"|). which to- 
gether with their Laplace transform determine the struc- 
ture of these operators, depend on the specific expression 
of the environmental state. It is therefore convenient to 
consider a specific choice, allowing for a more detailed 
evaluation. 



C. The vacuum case 

If the radiation field is in the vacuum state, the func- 
tions given in Eq. (TT3l) simplify considerably, since a(t) — > 
1, while (3(t) becomes a function of j(t) according to 
(3(t) — > \j(t)\ . The function ^(t) for the vacuum case is 
given by the expression 



G\t) 



— giAt/2 



COS 



~2 



A . 
%—— sin 



Vt x t 
~2 



,(32) 



where the superscript recalls that we have a single mode 
of the radiation field, while f^i = \J A 2 + 4g 2 . These re- 
sults for the vacuum case greatly simplify the expression 
of the obtained master equations, and inserted in Eq. (pH)) 
show that the time-convolutionless master equation off- 
resonance is always well defined. 

Indeed the time-convolutionless master equation for 
the vacuum case simply reads 



t^Vac 
'WCL 



(t) P 



-ilm 



-2 Re 



G J (*). 

GHt) 



[a + a-,p] 



a -pa j. 



(33) 



-{a+a-,p} 



corresponding as expected to a Lindblad structure with 
time dependent coefficients. Exploiting Eq. (I3"2l one has 



^TCh(t)P - 

9a 1 — cos (flit) 

19 A n! 



n ± t 

x [<r+(r-,p] 



fiit 



+ 2g 



2 sin(r2ii) 



cos 



Oii 



A 2 

n 2 



■ sin 



flit 



-i -l 



cr -P (7 + - 2 'f cr +' T -'^ 



, (34) 



where in particular one directly sees that the coefficient 
in front of the dissipative term at the r.h.s. of Eq. (|34|) pe- 
riodically takes on negative values, so that it describes a 
truly non-Markovian behavior and a Markovian approx- 
imation is not justified even for long times. 

The choice of the vacuum as bath state brings in 
important simplifications also for the expression of the 
Nakajima- Zwanzig memory kernel, whose Laplace trans- 
form reads 



!C$£(U)P: 



G\{u) 



G\{u) + G\{u) 

uz l (u] 



uz 1 (w) 



a -pa j, 



[a+a-,p] 



-{a+a-,p} 



GUu) 



z l {u) 



G R \u) + G]{u) i 



x [a z pa z ~ p] , (35) 



where 



z\t) = \G\t)\ 



(36) 



Taking into account the explicit result Eq. (|32|) one ob- 
tains after some calculations the following compact ex- 
pression for the memory kernel 



K^% c (r) p = -ig 2 sin (At) [a+a-, p] 



-NZ 



2g 2 cos ( y/A 2 + 2g 2 r 



(VA 2 + 2g 2 r 



<r-P<?+ - 2 { cr + (J -'/ ? } 



cos (At) [a z pa z - p] , 

(37) 



which is always well-defined even on-resonance. These 
results already allow for a few important remarks. We 
first notice that the different operator contributions in 
Lindblad form appearing in the various time local and 
integral kernels are multiplied by different time depen- 
dent functions. This lesson, already learnt in other mod- 
els [lj| [20, [29|, [30] , tells us that as a general rule non- 
Markovian dynamics cannot be obtained from Marko- 
vian expressions by simply taking the convolution with 



a single integral kernel, or multiplying by a single time 
dependent function. This fact has often been overlooked 
when seeking non-Markovian integrodifferential master 



equations (3U - I34J . sometimes leading to unphysical be- 
haviors. More than this, for the same model different 
sets of non-Markovian equations of motion can have dif- 
ferent operator structures, as it appears comparing e.g. 
the time-convolutionless and Nakajima-Zwanzig results 
for the vacuum Eq. (f33l) and Eq. (J37)) . This fact has al- 
ready been noticed in [20J , but the present analysis shows 
that this asymmetry depends on the choice of environ- 
mental state. For the present model it only appears in 
connection with the vacuum state. Indeed, while the dis- 
appearance of the term corresponding to excitation of the 
two- level system is obvious on physical grounds, when 
considering as bath state the vacuum, the vanishing of 
the coefficient in front of the dephasing term a z pa z — p 
is a peculiar feature of the time-convolutionless master 
equation. 



This model corresponds for a Lorentzian spectral density 
to the damped Jaynes-Cummings model. The time evo- 
lution map for this model, considering the special case of 
an environment in the vacuum state, has been obtained 

in [35] and can be expressed as: 



Pll (t) = pn(0)|G(t)| 2 , 

Pio(t) = P io(0)G(t), 

Poi(t) = poi (0) G* (t) , 

Poo(t) - Poo(0)+Pii(0)(l-|G(i)| : 



(40) 



where p(t) = <I> DVac (£)p (0), since we are considering the 
damped model with the bath in the vacuum state. The 
function G(t) is the solution of the equation 



IV. DAMPED TWO-LEVEL SYSTEM 

Exact master equations for bath in the vacuum 
state 



jG{t) = - J dti / (t - h) G fa) G (0) = 1, (41) 
with f(t) the two-point correlation function given by 



The technique used in Sect. IIHI to obtain the time- 
convolutionless and Nakajima-Zwanzig equation of mo- 
tions for a model whose evolution is known, by exploiting 
the representation of maps in terms of matrices, is appli- 
cable for a detailed study of simple models. We exploit 
it now for a variant of the model considered in Sect. [Til 
in which the environmental Hamiltonian is a collection 
of harmonic oscillators 



H T 



k 



and the interaction Hamiltonian is replaced by 
Hi = Y, (^ CT + ® bk + 9*k a - ® b l) ■ 



(38) 



(39) 



k 3 

= Y]\gk? e i{u °~ Uk){t ~ tx) , (42) 



corresponding to the Fourier transform of the spectral 

density. 



Starting from Eq. (|40|) and exploiting the same strat- 
egy used in Sect. Hill one immediately obtains for the ma- 
trix representation of the time-convolutionless generator 
the expression 



K^l c (t) = 



2 Re 



G(t) 
G(t) 



so that the master equation in operator form reads 



K-tcl (t)P 



-i Im 



-2 Re 



Re 
- Im 




6(d 



G(t) 
6(t) 

G(t) 







G(t) 
G(t) 

' G(t) 

G(t) 



Im 
Re 



[a+a- , p] 



CT-per-t 





G(t) 
G(t) 

6m 

G(t) 





1 



2 Re 



G{1) 



G(t) 



{<7 + cr_,p} 



(43) 



(44) 



which confirms the result obtained in [2J]- One can also determine the expression of the Nakajima-Zwanzig memory 
kernel, whose Laplace transform is given by 



jzDVac, 



U) 



( 



0_ 

Gr(u) 



uz(u)~- 1 



GrJu)+G t («) 

Gj(") 

G„ a («)+Gj 2 («) 



__0 

G H 2 ( U )+G, 2 (u) 
_ G«W_ 

Gr{u)+G! 2 {u) 





\ 



u£(u) — 1 



(45) 



where we have used the notation 

z{t) = |G(t)| 2 , 
leading to the master equation 



(46) 



IC^(u)p = +i- 



Gjiu) 



G R ( W ) + G 7 («) 



[(T+<T-,p\ 



1 — uz(u) 



z(u) 

1 — uz (u) 
z(u) 



< J -P (T + - -^{v+a-iP} 



2U- 



GJu) 



G R 2 (u) + G I 2 {u)}\ 

x \a z pa z - p] . (47) 



B. Role of correlation function 

We now show that this constructive approach leads to 
the same result as obtained by a suitable ansatz in [20l j . 
where the detailed connection with the perturbative ex- 
pansion via projection operator techniques has also been 
considered. We first have to recall the relation in Laplace 
transform between the functions G(t) and /(£). Indeed, 
for the case of a two-level system interacting with an en- 
vironment of oscillators in the vacuum state, all functions 
appearing in the non-Markovian equations of motion can 
be related to a single correlation function of the model, 
given by Eq. (|4"2"|). 

To fully exploit this fact in order to express the mem- 
ory kernel in the most compact and transparent way let 
us observe that according to Eq. (|4"Tj) one has 



- G*(u) 

f(u) = — j - U ' 



G(u) 



so that in view of the fact that for real 



(48) 



together with 



Re h(u) = Reh(u) 



Im h(u) = Imh(u), 



(49) 



(50) 



the functions 



fj{u) and - 



G/(u) 



g^,%) + g7 2 ( u ) 



(51) 



do coincide on the real axis, and therefore due to the 
identity principle on the common region of analyticity. 
A corresponding result holds for the functions 



/r(u) and 



G R (u) 



g^ 2 (u) + g7 2 (u) 



(52) 



so that the matrix representing the Nakajima-Zwanzig 
kernel can be rewritten in a more compact way as 



K 



DVac 
NZ 



(u) 



/ 



V 












\ 





-/*(«) 


~fi(u) 








//(«) 


-fn («) 





uz{u) — l 
z(u) 








uz(u) — l 
z(u) ) 



leading indeed to the integral kernel first obtained 



v 



Tz &c {r)P = -i//(r)[o- + <7_,p] 



(53) 
in(2fl| 

(54) 



-ki (t) 



< J -P< 7 + - 2 { a + a - 



>P} 



--{k 1 { T )-2f R {r)][cj z pa z 



P], 



where we have set 



fci(u) = 



1 — uz (u) 
z{u) 



(55) 



For the case of a single mode the correlation function 
considered in Eq. (|4"2"j) explicitly becomes 



!\t) = g 2 e 



2„iAt 



(56) 



where the superscript again stresses the fact that a single 
mode is considered. The solution of the integrodifferen- 
tial Eq. (|4"Tj) is then exactly given by the function G 1 (i) 
introduced in Eq. (f3"2"j). As it should be, the equations 
Eq. ((33} and Eq. J3J} are obtained from Eq. (JSJ and 
Eq. (|47j) under the replacement G{t) -> G x (t), which 
corresponds to the special choice of a single mode bath. 
At the same time it is clear that the explicit coefficients 
appearing in Eq. (|37[) can indeed be expressed using real 
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and imaginary part of Eq. (1561) . as well as the inverse 
Laplace transform of 



?1 1 - uz 1 (u) 



(57) 



2.9 2 



u 2 + A 2 + 2 ff 2 ' 
which is given by 

fcj(r) = 2 ff 2 cos(v / A 2 + V^) • (58) 



C. Additional term for thermal bath 

It is to be stressed, that the possibility to express 
all relevant functions appearing in time-convolutionless 
and Nakajima-Zwanzig master equations with reference 
to the single correlation function f(t) is a special feature 
of the two-level system coupled to the vacuum. Indeed 
while the results of Sect. IIIII for the vacuum are a special 
case of the model considered in Sect. IIV1 the more gen- 
eral situation of a thermal bath can be explicitly consid- 
ered for the case of the reduced dynamics of the Jaynes- 
Cummings model, showing that the expectation values 
of various operators as in Eq. (|13[) have to be specified in 
order to give the exact equations of motion. 

In Sect. IIIII we have seen that the time-convolutionless 
generator has a different operator structure with respect 
to the Nakajima-Zwanzig memory kernel only for the case 
of the vacuum, as it appears comparing Eq. (|33|l and 
Eq. (l37l) . while this is no more true for a thermal state. 
This strongly suggests that the asymmetry in the opera- 
tor structure of Eq. (|44p and Eq. (|54p is also due to this 
special choice of the bath state. To check this fact, in the 
absence of the exact time evolution map, one has to cal- 
culate the time-convolutionless master equation relying 
on the standard perturbative technique [2j, considering 
terms up to fourth order. The necessity to go up to the 
fourth perturbative order is immediately clear looking at 
the interaction Hamiltonian Eq. (|39|l . and observing that 
the dephasing term, as it appears from Eq. ([2"5]). involves 
a quadrilinear contribution in the raising and lowering 
operators <r+ and <7_. This task has been accomplished 
in Appendix lAl leading to the result 



^tcl(*)P 



ilm~f s (t) [cr + cr_,p] 



(59) 



+7+(*) 

+7_(t) 

1 



(T +P (T - - j {v-v+'P} 



<J-P(7- 



- 2 { a + a -'P} 



+ T7d(*) [°zPV z - p] , 

and the detailed expression of the various coefficients in 
terms of two- and four-points correlation functions of the 
system can be found in Eq. (JA22I) of Appendix [5] 



This result shows that indeed the disappearance of 
the dephasing term in the time-convolutionless master 
equation for the vacuum is a very special feature of this 
choice of bath. Once again the operator structure of non- 
Markovian equations of motion is strongly dependent on 
the details of both bath and interaction term. 



V. CONCLUSIONS 

We have obtained the exact time-convolutionless, time- 
local master equation and Nakajima-Zwanzig integrod- 
ifferential master equation for a two-level system cou- 
pled to a single bosonic mode, considering a generic bath 
state. This has been possible thanks to the knowledge 
of the exact dynamics, and corresponds to the result 
which can be obtained resumming all terms in the cor- 
responding perturbative techniques. The path followed 
here, to associate matrices to maps and to consider the 
direct relation between time evolution map and time lo- 
cal generator, as well as memory kernel, is however much 
more straightforward. The result shows that in a real- 
istic model each operator contribution in Lindblad form 
has its own time dependent function, responsible for the 
non-Markovian behavior, so that a simple multiplication 
or convolution of the Markovian result with a single phe- 
nomenologically guessed function will generally not work. 
Furthermore the operator structures of the equations of 
motion in the two cases can strongly differ, also depend- 
ing on the state of the bath. In particular it has been 
shown that a dephasing term, which is always present in 
the Nakajima-Zwanzig equations of motion, disappears 
for the time-convolutionless case for a bath in the vacuum 
state. This has been checked for the undamped two-level 
system thanks to the exact solution, and for the damped 
case by calculating the fourth order contribution to the 
time-convolutionless perturbation expansion. 

The pursued approach is quite straightforward, even 
though to obtain the exact expressions it relies on the 
knowledge of the full time evolution, which can obviously 
only be feasible for exceptional cases. The detailed anal- 
ysis of such cases however proves quite useful in under- 
standing the basic features of a non-Markovian descrip- 
tion, in particular it puts into evidence the strict rela- 
tionship between the different quantities which appear in 
the evolution equations, showing that phenomenological 
ansatz are in general not easily feasible. Furthermore, 
the connection to the microscopic perturbative deriva- 
tion techniques as outlined here, can help in determining 
the operator structure of the dynamical equations. This 
structure can sometimes be unveiled calculating the first 
perturbative contributions, thus restricting the number 
of ansatz necessary in order to consider a sound phe- 
nomenological model. 

The necessity to consider non-Markovian equations of 
motion is often due to strong coupling effects between 
system and reservoir. This in turn implies that a fac- 
torized initial state is not always the most appropriate 



initial condition, so that the study of the dynamics of cor- 
related initial states becomes of great significance. This 
long standing topic has recently seen important advance- 
ments [36I443J . and we plan to address it in future re- 
search work, looking in simple but realistic models, such 
as the one considered in this article, for general signatures 
of the effect of initial correlations. These results might 
help to consider, by means of some approximation, more 
general systems. 



J 
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Appendix A 



We here consider the contributions up to fourth order of the time-convolutionless projection operator technique for 
the damped two- level system considered in Sect. HVl We assume for the environment an equilibrium state pe, which 
commutes with the number operator, and we use the standard projection operator 



Tw = Tre (w) ® p E , 



(Al) 



where w denotes a state of system and environment. 

The time-convolutionless perturbative expansion can be obtained by repeated action of the projection operator 
considered above and of the superoperator 



C(t)w = -i[H T (t),w], 
where Hi(t) is the interaction picture expression corresponding to Eq. (|39|) . given by 

J3>(t) = a+{t)®B(t)+(j-{t)®B\t), 
with 



a ± (t) 



^±iu t 



and 



B(t) = Y, 9kb k e 



iuj k t 



(A2) 
(A3) 
(A4) 
(A5) 



The contributions to second- and fourth-order for the time-convolutionless generator then read [2] 



IC 2 (t) 



dtxVC{t)C{tx)V 



and 



/C 4 (i) 



dh r dt 2 [' dt 3 [PC(t)C(ti)C(t2Wa)V - VC(t)£(h)PC(t 2 )C{t 3 )V 

o Jo Jo 

-7>£(t)£(i2)P£(*i)£(*3)7> - VL{t)C{h)VC{tx)C{t 2 )V] 



(A6) 



(A7) 



respectively. Using these expressions one immediately obtains the operator form of the contributions to the time- 
convolutionless master equation for the reduced dynamics according to 



JC%l L (t)p = Tr E U dh 



ctftctydpstpE 



and 



/■*1 /"*2 

. |VL w > r ->'t \ I dtx I dt% I dt 3 

'o Jo Jo 



-(4) 



^tclWp = Tl 'J 



(A8) 



(A9) 



[P£(t)C(ti)£(t 2 )£(h)Vp ®p E - V£.(t)C(tiypC{t2)C(t z ypp ® pe 
-VC{t)L{t2)VC{tx)C{h)Vp ® P e- VC{t)C(t 3 )VC(ti)C(t2)Vp ® pe}} ■ 
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We recall that such high order contributions are needed in order to check the appearance, for an environmental state 
different from the vacuum, of the dephasing term a z pa z — p, which involves expressions with altogether four raising 
and lowering operators of the two-level system. 

The second order contribution Eq. (|A8|) can be expressed by means of the following two correlation functions: 



f(t-h) 



— Juo(t-ti) 



Tt E {B{t)B^(ti)p E } 



Y J \9k\ 2 e^~^ t - t ^{n k + l) E , 



(A10) 



which corresponds to Eq. (l4"2l if the bath is in the vacuum state, and 



— c — &>o(t— 1±) i 



g{t-ti) = e-"*v-*>Tt E {B\t)B(t{)pE) 



5> 



, 2 e -»(wo-w*)(t-ti). 



n k )i 



(All) 



which vanishes in the vacuum. In terms of these functions one has 

VC{ta)£(tp)Vp <8> p E = -[f(t a -tp)a+<j-p + f*(t a -tp)pa+a- 

+g{t a - tp)a-a + p + g*(t a ~tp)pa-a + 
-2Re/(t Q -tp)cr-p<j + - 2Rcg(t a -tp)a+pa. 



(A12) 



Pe- 



This result is sufficient to obtain the time-convolutionless master equation up to second order, indeed upon inserting 
Eq. (|A12[) in Eq. (|A8J) one obtains 



where we have set 



and 



KtclWp = -* [fi(t) + 0/ "(*)] W+<t-, p] 



(A13) 



+2f«(t) 



f(<) 



Jo 



cr-p(r+ - -{a+a-,p} 



dt x f{t-h) 



(A14) 



s(t) 



dtig(t- ti), 



(A15) 



denoting as usual real and imaginary parts with the subscripts R and / respectively. 

To consider the fourth order contribution one has to evaluate the four terms given in Eq. (|A9|) . The last three 
terms at the right hand side can be obtained applying twice the result Eq. (IA12[) . thus obtaining 



VC{t)C{t a )VC{t p )C{t 1 )Vp <8> pe 



[-4a + pa_{R,ef(t-t a )Reg(tp - i 7 ) +Reg(t - t a )Reg(tp -* 7 )} (A16) 

-4cr_ j ocr + {Re g(t - t a ) Re f(tp - i 7 ) + Re f(t - t a ) Re f(t p - £ 7 )} 

+<j+<j-pf(t - t a )f(tp - i 7 ) + pa + a-f*(t - t a )f*(t/3 - t y ) 

+<j-a + pg(t - t a )g(tp - i 7 ) + pa-(T + g*(t - t a )g*(tp - i 7 ) 

+a + a-pa + <7- {2 Re [f(t - t a )f*(tp - i 7 )] 

+4Rc 3 (i-t Q )Re/(i /3 -i 7 )} 

+(T_cr +/ 9(T_(T + {2Rc[g(t -t a )g*(tp -* 7 )] 

+4Ref(t-t a )Reg(tp-ty)} 

+<t-<t+p<t + <t- {f*(t - t a )g(tp - i 7 ) + g(t - t a )f*(tp - i 7 )} 

+cr + <T_ / 9cr_cr + {.g*(i - t a )f(tp - i 7 ) + /(£ - t a )g*(tp - i 7 )}] <g> ps, 



where the relations cr+ 



<r_ = have been repeatedly used, together with the assumption [p E , n k } = 0. 
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The first term at the right hand side of Eq. (JA9I) instead requires the introduction of a four-point correlation 
function, which is given by 



h(ta,tb,t c ,td) 



_ e iuo(t a — t b +t c -t d ) 



Tt e {B{t a )B\t b )B{t c )B\t d ) PE } , 



(A17) 



with B(t) as in Eq. (|A5|) . An explicit evaluation of VC{£)C(tx)C{t2)C(t 3 )Vp ® pe together with the repeated use 
of Eq. (|A16j) then leads to the desired result, which can be obtained with a straightforward though very lengthy 
calculation. The fourth order contribution reads 



£tcl(*)/> = * IM*) + «/(*) + «>/(*)] [<r+<r-,p] 



(A18) 



+t(t) 



+u(t) 






+ 4 [q(*) + *(*) + 2t)ii(t)] [a z pa, - p] , 



where in analogy to the notation of Eq. (|A14j) and Eq. (|A15j) we use the Fraktur character to denote the triple integral 
over time of the function with the corresponding Roman letter, for example 



pt pti pt 

p(t) = dh dte / 
Jo Jo Jo 



dt 3 p(t,ti,t 2 ,t 3 ). 



(A19) 



The functions determining the coefficients appearing in Eq. (|A18j) are given in terms of the above introduced two- 
and four-points correlation functions of the model according to the expressions 



p(t,t x ,t 2 ,t 3 ) = -^/(t-t Q )/(i /3 -i 7 ) + ^(i,i 1 ,i 2 ,i 3 ) 



(A20) 



afi'-y 

q(t,ti,t 2 ,t 3 ) = -2j2{ Re [f(t -t»)f*(tp -t-y)}+2Reg(t -t a )Ref(tp -t 7 ) -Reh(t a ,t,tp,t^)} 

r(t, ti,t 2 , t 3 ) = g{t - t 2 )g(h - t 3 ) + g(t - t 3 )g{ti - t 2 ) + f{t x - t)f{t 3 - t 2 ) - h(h,t, t 3 , t 2 ) 
s(t,t u t 2 ,t 3 ) = -253{Re[/(t-t a )/(fc r -t/9)]+2Re/(t-t a )Reg(t J 9-t r )-Refe(t,t a ,t r ,t /S )} 

a^7 

t(t,ti,t 2 ,t 3 ) = 2^{Rc[/(i-£ Q )/(i 7 -i /3 )]+Re[ 5 (i-i Q ) 5 (i /3 -i 7 )] 

+2Re/(£ — t a ) Reg(t/3 — £ 7 ) — Re/i(t,i a ,i 7 ,i / 3)} 

+2 {Re[/(ti - i)/(*3 - ta)] ~ Re[ 5 (i - h)g(t 2 - t 3 )} - Re h(ti,t, t 3 , t 2 )} 
u(t,ti,t 2 ,t 3 ) = 2^{Re/(t-t a )Re/(^-t 7 ) + 2Re 5 (t-t«)Re/(t^-t 7 )-Re/ l (t a ,t,t^,t 7 )} 

Ct/37 

-2Re/i(t,ti,t2,t 3 ) 
v(t,ti,t 2 ,t 3 ) = 2^2{f(t a -t)f{t r -t fi )-h{t a) t,t r) tfi)}, 

a/^7 

where the following summation convention has been used 

2j^(to:,t^,t 7 ) = i>{ti,t 2 ,t 3 ) + ip(t 2 ,ti,t 3 ) +ip(t 3 ,t 1 ,t 2 ). 

a/37 



(A21) 



Including terms up to fourth order one therefore has the expression Eq. (|59")l with time dependent coefficients given 
by the identifications 



nf a {t) = -f(t)-fl(t)+p(t) + t(t)+0(i) 

7+ (t) = 2f R (t)+t(t) 

7+(t) = 2g R {t)+u(t) 

ld {t) = q(t)+s{t) + 2t> R (t). 



(A22) 
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